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Abstract 

We introduce a minimal supersymmetric extension of a higgsless 
model for electroweak symmetry breaking in a warped extra dimen- 
sion. In contrast to the non supersymmetric version, our model nat- 
urally contains a candidate for cold dark matter. No KK-parity is re- 
quired, because its stability is guaranteed by an i?-parity. We discuss 
the regions in parameter space that are compatible with the observed 
dark matter content of our universe and are allowed by electroweak 
precision measurements as well as direct searches. 



1 Introduction 



The electroweak standard model (SM) has proven to be extremely successful 
in describing the hitherto collected data from particle physics experiments. 
However, we are still waiting for the forthcoming LHC experiments to find 

*e-mail: aknoch.el@ph.ysik . uni-wuerzburg . de 
^ e-mail: ohl@ph.ysik . uni-wuerzburg . de 



1 



the first direct experimental evidence for the detailed dynamics behind elec- 
troweak symmetry breaking (EWSB). In the minimal SM, EWSB is accom- 
plished by the condensation of a fundamental scalar Higgs field, but the 
corresponding particle has, so far, eluded detection at LEP and the Teva- 
tron. While the minimal SM with a rather light Higgs particle describes 
the electroweak precision data very well, any model with just a light fun- 
damental scalar is not natural when a much higher scale exists, such as the 
grand unification and quantum gravity scales. Therefore the minimal SM is 
expected to be just the low energy effective theory of an extended model in 
which new physics, that protects the mass of the Higgs particles, appears at 
the TeV-scale, e. g. supersymmetry (SUSY), little Higgs or extra dimensions. 

Instead of adding particles to the minimal SM in order to protect the 
naturalness of a light fundamental scalar, one can instead attempt to break 
the electroweak symmetry without invoking the condensation of fundamen- 
tal scalar fields. Indeed, such models were already proposed p] a few years 
after the introduction of the SM. In this class of models a new strong inter- 
action, called Technicolor (TC), effects EWSB in analogy to chiral symmetry 
breaking in QCD, but at the TEV scale. Unfortunately, after extending TC 
in order to produce fermion masses (ETC), the resulting models suffer from 
severe problems concerning flavor changing currents and electroweak preci- 
sion observable^- Recently a deeper understanding, based on the celebrated 
AdS/CFT duality, of such higgsless models has emerged which is inspired by 
EWSB in extra dimensions, particularly in warped spacetimes such as a slice 
of AdS5 [U [5]. Even though these approaches are nonrenormalizable and 
still suffer from some of the problems that plague ETC, they remain per- 
turbatively calculable up to energies of 5 ... 10 TeV [6] while simultaneously 
improving the fits to electroweak precision observables [TJ [H] ■ 

Combining observations of galaxy clusters, type la supernovae and the 
cosmic microwave background (CMB) from COBE and WMAP, there is now 
overwhelming evidence that nonbaryonic cold dark matter (CDM) consti- 
tutes roughly 20% of the energy density of the universe [9]. This CDM 
should predominantly consist of stable nonrelativistic, electrically neutral 
colorless particles (WIMPs). In light of this situation, it seems inevitable 
for any extension or alternative to the SM to eventually incorporate a phe- 
nomenologically acceptable CDM candidate. 

In SUSY models and in models with flat extra dimensions, there are nat- 
ural candidates for CDM that are stabilized by R- and KK-parity [TOl [TT] 
respectively. However, the warp factors and symmetry breaking boundary 
terms on the branes that are characteristic of higgsless models of EWSB 

1 Note, however, that some progress has been made in recent years 
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break the translational invariance in the extra dimensions. Therefore, the 
corresponding KK-parity is not a symmetry of higgsless models. As a re- 
sult, none of the Kaluza-Klein (KK) modes are protected against decay and 
they fail to provide a candidate for stable CDM. In the light of the above 
astrophysical observations, we are therefore compelled to extend the spec- 
trum of higgsless models in a different way if we want to retain them as 
viable alternatives to the minimal supersymmetric standard model (MSSM) 
and its cousins. One approach, which has been proposed recently [12], is 
to glue together two slices of AdS5 and to use an exchange symmetry as a 
KK-parity, protecting a CDM candidate from decay. Another approach [13] 
extends the spectrum and guarantees the stability of a CDM candidate by 
a Z 2 symmetry. 

While the naturalness of the TeV-scale does not require the introduction 
of SUSY in higgsless models, SUSY can nevertheless be expected to play an 
important role in any more fundamental theory. In particular, string theory 
requires SUSY at some scale. Indeed, string theory provided much of the 
inspiration for the revival of extra dimensional models, for higgsless models 
and for AdS/CFT. It is therefore natural to investigate supersymmetric ex- 
tensions of higgsless models. If these extensions admit an i?-parity, stable 
candidates for CDM are guaranteed to be part of the spectrum. Indeed, we 
propose to extend the above-mentioned higgsless models in a slice of AdS5 
with SUSY, providing a stable particle as a viable CDM candidate. 

In this paper, we introduce the following framework: assume that the 
sector providing EWSB exhibits a global SUSY whereas the SU (2)lXU (l)y 
gauge theory to which it is coupled does not. Our approach extends the par- 
ticle spectrum and its interactions in a well-defined way requiring a minimum 
of new assumptions. We examine possible scenarios which differ in how the 
fields on the UV brane are coupled to the bulk, calculate the relic density 
and show that the observed dark matter density Qh 2 « 0.09 . . . 0.11 [9j can 
be achieved straightforwardly with reasonable gaugino masses. 

The paper is organized as follows: in section [2] we fix our notations, by 
giving a concise review of SUSY and the SUSY spectrum in a warped back- 
ground. In section [3] we outline the supersymmetric version of the higgsless 
models on which our construction is based [UEHHIH]. In section H] we survey 
the parameter space of our proposed model, identifying the CDM candidate 
and estimating its relic density. We find that it is possible to obtain realistic 
relic densities of a neutralino LSP without much fine tuning. Section [5] sum- 
marizes our conclusions. Further technical details, regarding the interplay of 
SUSY, KK-decomposition and boundary conditions (BCs) can be found in 
appendix |A] 



3 



2 Global Supersymmetry in a Slice of AdS5 



In this work we assume that the backreaction of the matter and gauge de- 
grees of freedom — including superpartners and KK-towers — on the warped 
geometry can be neglected. Therefore, we only need to implement the global 
SUSY transformations that are compatible with the isometries of the warped 
background geometry and not the full supergravity algebra. 

In the case of flat extra dimensions, the SUSY algebra in Minkowski space 
can be generalized straightforwardly to five dimensions by using the corre- 
sponding Clifford algebra and promoting the parameters to four component 
spinors. The commutator of two transformations then still reads 

[5 n ,8 ( ] = -2(irf i t-tT u Ti)Pu, (1) 

where the gamma matrices are now defined such that 

{ 1 M , 1 N } = 2 V MN , M,iV = 0...3,5. (2) 

Expressing the 5D iV=l-SUSY algebra in terms of four component genera- 
tors, 

{Qi,Qj} = -HJPm = -2 7 g^ - HA , (3) 

a comparison with the 4D iV=2-SUSY algebra reveals that iP$ plays the 
role of a central charge in the 4D picture [14] . Since translations leave the 
metric i]mn invariant, we are dealing with a global spacetime symmetry. 

Moving on to a warped background, we retain the approximation of ne- 
glecting the backreaction, as stated above. Following the approach of [15] for 
the treatment of a curved background, we define global SUSY transforma- 
tions by demanding that they leave the metric g invariant. In other words, we 
demand that these global SUSY transformations close into a Killing vector 
field v of the background metric. In particular 

[6t,6 v ] = v N P N (4) 

and the Killing condition for v reads 

v M d M g AB + g A Md B v M + g B Md A v M = , (5) 

while the gamma matrices satisfy now 

{7^1 = 2/™. (6) 

This gives us a condition for the spinor valued parameters of the SUSY 
transformations, and solutions (£, rj) of (TjO) and (jHJ) are called Killing spinors. 
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We use a "mostly — " metric convention. The warped background metric 
we are dealing with in this work is that of Randall-Sundrum (RS) Type I [16J, 
which, in proper distance coordinates, reads g^ v = e~ 2Rfcj/ ?7^, #55 = — R 2 , 
and consequently y/g = Re~ 4Rky . Here rj^ denotes the Minkowski metric, 
R the radius, k the RS curvature, and y G [0, tt]. As customary, the 4D 
spacetime at y = (y — tt) will be referred to as the UV (IR) brane. We 
work in the interval picture with BCs, noting that this space could also be 
interpreted as an orbifold. There are the usual 4D Poincare symmetries and 
an additional scaling symmetry x M — > (1 + 5)x M which is broken only by 
the presence of the branes. Working out the Killing condition (pH) and (jSJ) in 
this background, one ends up with a set of SUSY parameters which generate 
SUSY transformations that close into the remaining symmetries, namely, 
using 2-spinor notation, 



£(x,y) = e-^ , (7) 



to 

So: 



where the space-time dependence remains confined to the warp factor. This 
relation fixes the KK wavefunctions of the superpartners. Since (J7|) is param- 
eterized by a single Weyl-spinor, there can be at most one 4D supersymmetry 
left after integrating out the extra dimension. Nevertheless, we will see in the 
following sections, that the spectrum of the massive KK modes will formally 
be that of 4D iV=2-SUSY. 

Counting the degrees of freedom reveals that a massless 5D vector boson 
can not be combined with 5D spinors to form a SUSY multiplet. The smallest 
multiplet that contains both is actually a dimensionally reduced 6D vector 
multiplet which consists of a 4D vector multiplet and a chiral multiplet p2] . 

On the other hand, the 5D hypermultiplet is constructed from a Dirac 
spinor ^ = (ijj a ,ip ca ) T and two complex scalars as its superpartners. This 
makes it equivalent to a combination of one 4D chiral and one 4D antichiral 
multiplet. These relations have been extended to a full 4D iV=l-superfield 
formulation of these multiplets [T7J [18], where the 4D kinetic terms arise 
from the superfield itself as usual while the 5D dynamics are put in explicitly. 
Following [T5] , the superfields take the form 

V a = -6a m 6A a m - ite\ a x + i6 2 6\l + -6666 D a (8a) 

X a = ^^ a + iA t) + -^fa m 6d m {Y> a + iAf) - -l=^MD(S a + iA a 5 ) 

+ y/26\ a 2 - -^66d m X a 2 a m 6 + 66C a (8b) 
V2 
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for the gauge multiplet in Wess-Zumino gauge with the adjoint index a and 

H = h + i9a m 9d m h - -eeWUh + y/28ip - -^ r 99d rn ipa m 9 + 99F (9) 

for the chiral superfields H and H c in the hypermultiplet. In this notation, 
the action of the gauge multiplet for the nonabelian case is given by 

S g [V, X ] = j d 5 x J d 2 9 ^ tx[W a W a ] + h.c. 

/r „—2Rky 
d 5 x j d 2 9 ti[(V2d y + x)e- v (-V2d y + xV + d y e~ v d y e v ) , 

9 (10) 

where V = V a T a and x = x a T a are the Lie- Algebra valued superfields with 
normalization tr(T a T b ) = N5 ab and g/2 is the 5D coupling constant. The 
bulk action of the hypermultiplet coupled to the gauge fields reads 



S h [H, H c , V, X ]= J d 5 x J d 2 9 Re~ 2Rky 



He~ v H + H c e v H c 



+ Jd 5 x J d 2 9e- 3Rky H c 



0„ 



V2 



X 



c)Rk 



H + h.c. (11) 



Note that d y , x an d V are dimensionless whereas H and g have mass di- 
mension 3/2 and —1/2, respectively. The dimensionless quantity c is the 5D 
mass of the multiplet in units of the RS curvature k. After the field redefini- 
tions described in appendix [B], we obtain canonical kinetic terms. With our 
conventions, the gauge couplings in this action correspond to the covariant 
derivatives 



D M V = {V M -iA a M T a )^ 
(D M $) a = V M $ a + f ahc A b M ^ c . 



(12) 
(13) 



Using these ingredients, we can now go on to construct the supersymmetric 
model. 
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3 The Model 



3.1 Higgsless Electroweak Symmetry Breaking 

Let us first give a short description of the Higgsless models our construction 
is based upon [5j [19] . They feature a left-right symmetric gauge groupEl 

G = SU(3) C x SU(2) L x SU(2) R x U(l) x (14) 

where X will be the (B — L)/2 quantum number which is 1/6 for quarks and 
— 1/2 for leptons. The corresponding coupling constants are g 5s , Ql = 9r = 95 
and (75. The symmetry breaking by BCs is designed such that 

(SU(3) c x SU(2) L xU(l) Y ontheUVbrane 

\SU(3)c x SU{2) D x U{l) x on the IR brane "' >] 

leaving only an overall SU(3)c x U(1)em unbroken. The diagonal subgroup 
SU{2) D generated by Tg = T£ + T£ acts as custodial symmetry. The BCs for 
the gauge fields (and later on those of their scalar and fermionic superpart- 
ners) must be compatible with this breaking scenario and with the vanishing 
of the variation of the boundary action. This is discussed in appendix [B] 
together with the resulting KK decomposition. By virtue of these BCs, the 
W ± is a mixture of the SU (2) L and SU {2) R gauge bosons (localized in the 
SU{2) L field up to one permille in order to ensure the observed V — A cou- 
pling), the Z and 7 are a mixture of the SU(2) L , SU(2) R and £7(1) gauge 
bosons. The masses of W ± and Z are given approximately b}{§ 

U e -2Rkir 

m ™ " (TT^ (16a) 

2 (g 5 /g 5 ) 2 + k + 2 2 

to/w + 1 mw ' (16b) 

(where the brane kinetic term ( 1531) contributes the ^-dependence) so the 
radius R is determined by the RS curvature k and the W mass. 

The leptons and quarks are implemented as in [7] . There are two doublets 
transforming under SU(2)l and SU(2) R respectively for each SM fermion 

^ L = (r L ,W,<W c ) T (17a) 



2 The general case 7^ has been studied in the literature and it turned out not to 
be an effective means to improve precision fits and perturbative unitarity. We therefore 
assume the gauge action to be LR symmetric in the bulk for simplicity. 

3 If the ratio 35/^5 is determined from the particle masses, the deviations from the SM 
are shifted to the couplings. To make the tree level corrections "oblique", this quantity 
should be defined by fixing the gauge couplings to matter first [?]■ 
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y R = (r R ,r R c ,« c ) T . (17b) 

The two doublets get 5D Dirac masses denoted by cl and cr respectively 
(cf. ( Hip ) which are allowed in the bulk where the theory is vectorlike. In 
the limit of massless fermion zero modes, we impose BCs such that ipL and 
ip c R have a zero mode for which ifjR and ipl vanish. The mass of the zero 
mode is then lifted by a SU (2) ^-invariant Dirac mass fi on the IR brane, 
resulting in modified effective BCs which mix and ^r. So far the treat- 
ment is the same for the quarks and leptons. They differ however in the 
way how the doublets are split on the UV brane where the theory is not 
invariant under SU(2)r,. Quarks are given a UV brane kinetic term with 
parameter p, whereas the SU{2) R neutrinos can receive a large UV localized 
Majorana mass p r which leads to a seesaw-like mechanism. The mass of 
each fermion is determined by the UV splitting parameter, the IR mass and 
the 5D "bulk" masses of the doublets, cl and cr, which control the localiza- 
tion of the fermion zero modes. The resulting masses are approximately (for 

C L > 1/2, Cr < -1/2) 

y/(l-2c L ) (l + 2c R ) (e^Y^- 1 p 



m,, 



/ l-(l + 2c /? ) hp* 

CL—CR — l 



1) (e' 2 ^) 



18a 



(18b) 



3.2 The Supersymmetric Extension 

The structure of the 5D multiplets outlined above implies that we must 
promote all fields to superfields separately 



if, 



V 
X 
H 
H c . 



Now we can write the supersymmetric BCs for the gauge multiplets in a 
form on the IR brane (i. e. y = it) 



(19a) 

(19b) 
(19c) 

short 
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and on the UV brane (i. e. y 
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Neutral Charged Colored 

=1 S=l/2 S=0 S=l S=l/2 S=0 S=l S=l/2 S=0 



Figure V. The particle spectrum with unbroken SUSY for a typical set of 
parameters Rkn ~ 37.5 (k ~ Api). The S = 1/2 bars each represent 
one Majorana particle in the neutral column and one Dirac particle in the 
charged column, while the filled/unfilled 5 = 1 bars distinguish between 
massive/massless vectors with three/two physical degrees of freedom 

d y V L (0) = U m2 (0) = d y V c (0) = (20d) 
X L (0) = d y e- 2Rk y X R12 (0) = X C (0) = . (20e) 

These BCs are satisfied by solutions to the Euler-Lagrange equations derived 
from the free part of (1101) . The BCs of the canonically normalized fields are 
obtained after performing the redefinitions (H71) in appendix [B] As it has 
already been noted by the authors of [HI ED], the tree level spectrum re- 
mains highly degenerate even though one of the supersymmetries has been 
broken by the warped background (shown for our case in Fig. [T]). Similarly, 
the generation of matter fermion masses on the IR brane is made supersym- 
metric by giving SU(2)o invariant localized masses to the entire multiplet 
5(tt — y){H1Hn + HfjHi). Like for the case with fermions only, there are var- 
ious prescriptions to convert this into effective BCs, and the hypermultiplet 
becomes degenerate for unbroken SUSY. 

Of course, similarly to the MSSM, such a spectrum is not acceptable be- 
cause massless neutralinos, gluinos and light charginos below 95 GeV have 
already been ruled out by experimental searches. Also, the gauge scalars 
are even under R parity and could be produced below the pair production 
threshold. We will therefore investigate the consequences of having no super- 
symmetry at all on the UV brane. Starting from this premise, it is natural 
to remove all scalars from the UV brane by imposing the BCs 

S L (0) = £*(0) = S x (0) = S c (0) = V L (0) = /4 C (0) = . (21) 
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This pushes the SU(3) and U(l)x gauge scalars up to a masa^ 



m s c = m s x = z ke~ Rkw ~ z Q ^/Rkn{l + k) m w ~ 1.2 TeV (22) 

where zq ~ 2.41 is the first zero of the Bessel function Jq(z). The sfermions 
lie in a similar mass range, depending on the localization parameter c of the 
multiplet. For a massless fermion we can approximate the tree level sfermion 

masses 

m; = z ke~ Rkn , (23) 

where, if cl ^ 1/2, cr ^ —1/2, zq is the first positive root of J CL -i/2 for Jl, 
J\c r +i/2\ for fn, J\ CL -i/2\ for /£ and J- CR -\/2 for /£. The situation is different 
for the scalars from gauge groups that are broken on the IR brane. These 
scalars receive smaller tree level masses 



m E o = m s+ ~ J -j^e- Rkw = y/2(l + k) m w , (24) 

which is very interesting from a phenomenological point of view. However, 
their tree level coupling to fermions is of the form Yiipij) and therefore sup- 
pressed with the fermion mass, vanishing altogether in the massless fermion 
limit where it is forbidden by the chiral symmetry. For a large range of pa- 
rameters, the coupling to leptons and quarks is similar to the corresponding 
SM Higgs coupling. Note however that since it is not the vacuum expectation 
value of this scalar triplet that breaks electroweak symmetry, this similarity 
does not extend to the coupling to gauge bosons. Consequently, the processes 
corresponding to Higgsstrahlung and vector boson fusion in the SM are ab- 
sent in our model at tree level. This provides an experimental signature for 
distinguishing our model from the MSSM at the upcoming LHC experiments. 

Now let us turn to the gauginos. We split the gauginos and the gauge 
bosons with the localized kinetic term ( l53li 5 l. We will study the following 
two sets of BCs which project out all massless modes, while maintaining 
the SU{2) L x U(l)y invariance (but not the full supergauge invariance) and 
the stationarity of the action on the UV brane. In the first scenario 

\% = Af = Af = , (25a) 

4 This is accurate if further localized kinetic terms are absent. The values are given 
here for k = Ap\. 

5 If one does not mind tachyonic solutions above the cutoff of the effective theory at 
~ 10 TeV, it is also possible to raise the chargino mass sufficiently with a localized kinetic 
term, but this possibility will not be investigated further in this work. 
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Figure 2: The tree level contribution to the S parameter for a massless probe 
fermion, Rkir ~ 37.5, k ~ Api , cr = — cl and different values of k. (k = 
0, 0.2, 0.4, 0.6, 0.8 from highest to lowest) 



while in the second 



Af = Af 2 = 



cos(#)Af + sin(0)Af 3 = cos(#)Af - sin(#)A* = . 
In both cases, the lightest charginos will receive a tree level mass 

Vl + nm w , 



(25b) 



(26) 



which requires k > 0.4 to escape the current detection bounds. The resulting 
mass of the lightest neutralino mode will then be 



m x o = m x -i 



in the first scenario (125aj) and 



m x° — 



cos(#) m y 



(27a) 



(27b) 



in the second (I25bl) . 

If we isolate the oblique corrections that arise at tree level analogous to 
[211 [7] , there are two parameters that have a strong influence on the param- 
eters S, T and U in our scenarios, namely k and the fermion localization 
parameters c. Using the method of these authors, the resulting value for S 
is given in fig. [21 We find a combination of the results in [7] and [21] • The 
kinetic term used to split the electroweak gauge multiplets moves S in the 
right direction and can allow for delocalized or slightly UV localized fermions. 
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3.3 Supergravity 

In some supersymmetric extensions of the SM, the gravitino poses a prob- 
lem for cosmological observables, either because it is overabundant or so 
long-lived that its decay products spoil nucleosynthesis. So far we have not 
included supergravity in the bulk, and we want to give an argument that 
the inclusion of KK gravitons and gravitinos in the bulk is possible with- 
out loosing the desired properties of the model. A holographic interpretation 
along the AdS / CFT correspondence where the energy momentum tensor and 
the supersymmetric currents are present on the CFT side also warrants the 
inclusion of supersymmetric gravity multiplets. Similarly to BC gauge sym- 
metry breaking we adopt the philosophy that the simple picture with BCs is 
a description mimicking a full theory with symmetry breaking in a slice of 
AdS5 [15j [22] • By setting the appropriate BCs we obtain a heavy short-lived 
gravitino. 

Solving and KK expanding the (free) Rarita-Schwinger equations with 
"twisted" IR BCs ^2(0) = ^li 71 ) = 0, we obtain a very light gravitino 
m — ^/8ke~ 2Rk7T [20] . However, twisting the BC on the UV brane yields the 
condition 

T CHL\ T ( p Rkitm\ 

Y 2 (f) Y^f) rn n -Z n Ke , 

where the positive roots of J\. From zq ~ 3.83 we find a heavy 

gravitino solution at the scale of the lightest graviton KK mode (0(2 TeV) 
for k = Api). Like the heavy gravitons, this gravitino does not have Planck 
suppressed interactions but rather a coupling < TeV" 1 depending on the 
localization parameters. To be more precise, let us calculate the coupling 
scales relevant for neutralino annihilation (and gravitino decay). The general 
solutions for the KK wave functions are (cf., e. g., [20, 1221 [23] ) 



GT\y) = e 2Rky [aJ 2 {m n e Rky /k)+(3Y 2 {m n e Rky /k)\ (28) 

Rm r 



p—Rky 

& 2 n \y) = e 2Rk y [aJ 1 (m n e^/k)+0Y 1 (m n J a "/k)] = <9,G< n) , (29) 



' n 

with the canonical normalization condition 



With these conventions, the coupling strength to a vector and a gaugino is 
(GifvU = ^{f k J ^ ,mky Gifvf Xi ■ (30) 
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The normalization is chosen such that a gravitino zero mode with unbroken 
SUSY would yield (G { ^f v f Xl ) = A" 1 , G< 0) = 0. For « = 0.4 and k = A P1 we 
obtained 





a (500 TeV)- 1 




(25 TeV) -1 


(31a) 




a (lOTeV)" 1 




(80 TeV)" 1 


(31b) 


(Gi/!°/$> p 


a (lTeV)" 1 


(G 2 /f /|) « 


(3 TeV)" 1 . 


(31c) 



Considering the usual holographic picture, this result is not surprising, be- 
cause the heavy gravitino is interpreted as a bound state with TeV suppressed 
interactions to other bound states f)31cl) . coupling less strongly to lighter, 
mostly elementary particles (131 al) and (131bl) . With the result in f)31al) for the 
SM couplings we can neglect neutralino annihilation with the gravitino in the 
t-channel in the following discussion. Still, such a gravitino couples strongly 
enough to essentially vanish immediately for a temperature of T 2 TeV, 
long before nucleosynthesis. 



4 Dark Matter Phenomenology 

Our aim in this section is to estimate the relic density of the neutralino 
LSP we found in the scenarios discussed above. The analysis of thermal 
WIMP production and relic densities generally requires to numerically solve 
the Boltzmann equation for the scattering processes involved, but there are 
approximate semi-analytic solutions adequate for WIMPs in SM extensions. 
In our calculation of the neutralino freezeout temperature and relic density 
we follow [UJ ITU] . In the case without coannihilationa^l we need the number 
of effectively massless degrees of freedom g*, the LSP mass and the total 
cross section x^x'j ~^ A expanded in the relative velocity v 

valf(v) = a ij + b ij v + 0(v 2 ), (32) 

which we obtain by evaluating the graphs in Fig. El The cross section (I3"2"j) 
can be reexpressed in terms of the temperature T 

(vaif)(T)^a t3 +b z3 ^. (33) 



6 Even though the neutralinos may be split by radiative corrections we treat them as 
one particle with four d.o.f. rather than having coannihilation between them. 
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Figure 3: Tree level contributions to neutralino annihilation. 

The relative freezeout temperature 1/xf — Tf/m x o can be determined itera- 
tively, and is typically at around Xf ~ 15 . . . 25. The final expression for the 
relic density is then 

_, 2 , , 1.04 x K^GeV" 1 x f 1 , g .. 
Api y/g*a + 3b/Xf 



The effective number of degrees of freedom is taken to be g* = 86|, to which 
fermions contribute with a weight of 7/8 from the Fermi distribution. To 
give an estimate of the individual contributions XiX'j ~^ X n , we will add the 
inverse densities 

l/tth 2 (x f ) = tthKxf)- 1 oc a «+ 3 b n/ x f ( 35 ) 



for a realistic fixed freezeout temperature, which depends only weakly on 
the cross section. The value currently favored by WMAP data is 1/Qh 2 « 
9. ..11 [S]. 

4.1 Neutralino Annihilation Channels 

Depending on the neutralino mass there are several annihilation channels 
contributing at tree level (Fig. [3]). The channels XiX'j ~ > ^> p v an d XiXj ~~ * 
dd, uu, ss, cc, bb — > jets are always open, whereas Xi X°j ~> W + W~ is open 
for all parameters in scenario (I25ap but only if cos^i/l + k > 1 in sce- 
nario ( I25bl) . We first rewrite the neutralinos (charginos) as two Majorana 
(Dirac) fermions. To achieve this, we express the 4D action in terms of the 
KK coefficient fields \\' 2 (x), X^' (x), A° (x) and A° (x) and group 
those into 4 spinors (the KK indices are implicit) 



x^) = ^(a? + a°,a? + a°) 
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Zn x 1 




z,i x' 




Figure 4: Potential tree level contributions to neutralino annihilation from 
the gravitinos and massive gravitons 



X&0 = ^(A?-A§,Ag-A? 
X t{x)^^{x\-iX\X-ixt) 

With these definitions, the full 4D chargino-matter-smatter interaction La- 
grangian takes the form 



C = h 5 ViiTM + ixl) + hf VtiTUxa ~ ixl) + h.c 



(36) 



The vertices are given in terms of the projectors, overlap integrals and quan- 
tum numbers by 



+ ~g 5 X, tl [p-{fx-h Li h L] ) -P + (f^U h fh Lj ) 
+ g 5 X i:j P~ (f x x f i>m f hR . ) - P + (f x x f^ m f h 



Rj 1 



-Rjl 



r-- = <,:Jf [p-(h?UJ h% ) + P + {h?hiJ hl] ) 
+ 9^ R [P~{fxpU m fh Rj ) + P + (fx^fr R Jh% J ) 

+ g 5 X ZJ [p-(f^U u h h ) + P + (f xX UcJ hlj ) 
+ g 5 X i:j P~ (f x x f lPm f h c R , ) + P + (f x x f Vm f h c R . ) 
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In this expression, the brackets (fxUfh) = J dy^/g fx{v)U{v)Sh{y) stand for 
the 5D overlaps which give us the coupling strengths of a sfermion to the 
corresponding matter fermion and a neutralino. The structure is such that 
the "left handed" sfermion couples to gauginos and matter fermions of the 
same handedness, the "right handed" one to gauginos and matter fermions 
of opposite handedness. The charged current interaction Lagrangian is 

£ = [^+ C {X° a ~ iX° b ) + xfK (X° a + »X»] + h-c- (37) 

and the corresponding vertex expressions are 

r a = 95l n P~((fw L ±fx^±fx^) + (fw R ±fx R± f\? 3 )) 

+ p+ {(fw L ±f\[±f\Y A ) + (fw R ±fx R± fx R3 )) 

= r n a (p + ~ p-) , 

with the 5D overlap integrals (fwfxjx,) = Jdy^/g e Rky f w (y)f Xl (y)fx J (y) 
which now have an additional factor e Rhy from the inverse vielbein contained 
in 7 M . The effective coupling constants for the charged current can be approx- 
imated analytically to leading order using Rkn 3> 1 and my/ 1 m x -C ke~ Rkn 
as 

3 1 

(fw L ±fx^fx^) — -gC (fw R ±fx R± fx R3 ) — -g C ' ( 38 ) 

with 

(fw^fx^fxf 3 ) - ~ C (fw R ±f x R± fx R3 ) - ~ A8RklT C ^ 39a ^ 
in scenario (I25al) and 

(f W L±f xL± f xL3 ) ~ J-^pic (f wR± f xR± f xR3 ) ~ - 7|c0 " (g)l C (39b) 
\jw j x r a, / 2ARkn XJW Jx i Jx * 1 v ; 

in scenario ( 125bl) . respectively, where 

C = -. 1 . (40) 

y/(K+l)R7T 

They are accurate to about 5% for Rkn > 30 (cf. also Fig. [3] below). Note 
that the absolute size of the coupling is approximately independent of the 
free parameters, 

n v/ 4™; ... 
9 5 C ~ . , 41 

and the leading contributions are also independent of the neutralino mass. 
For the sfermions with mj > ke~ Rkn the same expansion is not possible. 
The Feynman rules were implemented and the cross sections evaluated using 
FeynArts/FormCalc [2i 
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4.2 Results 



4.2.1 Annihilation into Fermions 

Due to the large KK mass of the sfermions around — cr, cl ~ 0.2 . . . 0.7 of 
mj > 2ke~ Rk7T , the cross sections x°X° f f wm De suppressed relative to 
the annihilation to W pairs if m x o > m^. We find that the couplings and 
cross section depend strongly on the localization of the fermion multiplets, 
but that a rather extreme localization towards the TEV brane is necessary to 
make a significant contribution to the total annihilation cross section. This 
region of parameter space is ruled out due to FCNCs (unless an additional 
flavor symmetry is introduced to suppress them) and proton decay from 
higher dimensional operators. Therefore we conclude that in this scenario 
the tree level annihilation of neutralinos into fermions is negligible. For 
example, if one carries out the calculation for k = Api, k = .4, Xf = 23 
the contribution of the charm quark stays below 1/Qh 2 < 0.3 for \cl,r\ > 
0.25 and even 1/Qh 2 < 0.02 for \c^p\ > 0.5. The electron which must 
be at cl ~ 0.5 to get a realistic S parameter, contributes 1/Qh 2 < 0.01. 
One uncertainty comes from the precise implementation of the third quark 
generation. The straightforward generation of mb and m t with boundary 
terms leads to problems with the Zbb coupling^. Depending on how the b 
quark is realized, its coupling to the neutralino could be somewhat enhanced, 
but we can still expect that the channel XiX^ ~~ * bb does not play a major 
role for dark matter annihilation. 

4.2.2 Annihilation into W Pairs 

If m x o > mw, the neutralino LSP can annihilate into two Ws at tree level, 
and the interaction is essentially of weak strength (jSJ). In scenario (I25aj) the 
annihilation cross section turns out to be too large, leading to a small relic 
density which would require another mechanism for Dark Matter production. 
In scenario ( I25bl) the largest coupling constants are strongly suppressed, and 
the relic density is in a realistic range (Fig. [5]). A neutralino LSP with 
""M ^ 86 GeV is favored. 

5 Conclusions 

We have constructed a minimal supersymmetric extension of warped higgsless 
models. The absence of SUSY on the UV brane leads to a tree level spectrum 
in which the lightest neutralinos are naturally the lightest superpartners and 

7 For a proposed solution see [8] . 
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84 86 88 90 92 

mx°/GeV 



Figure 5: The m x o-dependence of contributions to 1/Qh 2 for annihilation 
into W pairs for Rkir ~ 37.5, k = Api. The parameter is (from topmost) 
k = 0.4, 0.6, 0.8, 1.0, corresponding to m x + = 95, 102, 108, 114 GeV. For this 
process we use a = 1/128. The dashed lines are obtained using the approxi- 
mations in (I38 j) -( I4TI) . The shaded area represents the region allowed by the 
WMAP three year data. 



the lightest charginos the NLSPs. Experimental constraints force us to shift 
the chargino mass above 95 GeV, with m x o = . . . m x + depending on the 
mixing angle. All color charged superpartners and the sfermions are at the 
KK scale < O(TeV). 

The lightest electroweak gauge scalars S 1 * 1 , S° have tree level masses 
at V2m x +, but are expected to receive significant radiative corrections. Their 
Higgs-like masses and couplings to matter make them an interesting object 
for searches at the LHC, where they could be distinguished from the Higgs 
by the absence of Higgsstrahlung and vector boson fusion processes. These 
issues will be treated in an upcoming publication, along with a more detailed 
analysis of LEP precision observables. 

Another aspect to be investigated is tree level unitarity of the model which 
generally requires vector boson KK modes lighter than the ones found in our 
minimal scenario, a feature which can for example be achieved by including 
additional boundary terms. 

For m x o > mw, the dominant channel for LSP annihilation is into W 
pairs, while the annihilation into SM fermions does not contribute signifi- 
cantly unless they are highly localized towards the TeV brane. If the left- 
handed neutralino is set up to be localized mostly in the SU(2) R triplet, 
the LSP relic density is in agreement with current observations for a nat- 
ural range of parameters where m^y < m x o < m x +, making SUSY in the 
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EWSB sector a promising source for the observed nonbaryonic dark matter 
if Higgsless symmetry breaking turns out to be realized in nature. 
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A Conventions 

We use the flat metric convention rj = diag(+, — , — , — , — ) and the corre- 
sponding Dirac matrices 

where o~° = <t° = — 1 and —a 1 = a % are the Pauli matrices. We define the 
projectors 

P+ = iH 7 § + l) P- = 1h 7 5 -1) (43) 

on the 4D "left handed" and "right handed" component respectively. For 
Dirac spinors \1/ T = (i] a ,x a ) we define \I/ = = (x a ,Va)- 

Two coordinate systems are commonly used. The "proper distance" coor- 
dinates have x M = (x^, y) and y G [0, ir}. They are related to the "conformal" 
coordinates with x M = (x^^z), z e [A; -1 , A^ 1 ] through 

z = k~ 1 e Kky Air = ke~ Rk7r . (44) 

The Dirac matrices read in proper distance coordinates 

7 ^ = -^e^(_° m ^) (45a) 



^-r[o -i ) W 
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7„ = V,»e- 2Rky Y (45c) 
7 5 = -i?V (45d) 

(45e) 



and in conformal coordinates 



/ cr m \ 

-f = -S&kz[£ m a Q J (46a) 

7 5 = ^(q ^) (46b) 

7m = ^^7* (46c) 

75 = -i7 5 - (46d) 



B Boundary Conditions and KK Decomposi- 
tion 

Boundary Conditions 

The kinetic terms that we obtain straight from the superfield action and the 
BCs take the canonical form after the field redefinitions 

(47) 

A 5 — > gRA 5 A 2 — ► -igRe-2 Rky \ 2 £ — ► 

The standard coupling constants are recovered after taking g>5 = g/2. Before 
we specify the KK decomposition for the fields, we employ an type gauge 
fixing term 

3 gf = -J d 5 x — l V mn d m A n - t^-{dy ~ 2Rk)A 5 ) . (48) 

The BCs in ( 1201) now give the KK modes A^ an unphysical mass \f^m n if 
the corresponding massive gauge boson mode has a mass m n . This can be 
seen as follows: With (j4"8|) and (jTUl) the KK equations of motion are 

m 2 A( n) = e -2Rky ( _ d 2 /R 2 + 2^/^)4") 

Taking the derivative of the first equation, we find 

m 2 n d y Af = e~ 2Rk y(-d 2 y /R 2 + Akd y /R - Ak 2 )d y A^ . (49) 

20 



Since d y A^ satisfies the same equation of motion as A 5 , the correct choice of 
BCs immediately follows: 

dyA^g = => A 5 \ 9 = 

m ^ (50) 

A M | 9 = =4°(^-2i^L4 5 |s = 0. 

Now, there is a mode of A 5 for every massive mode of A^. A similar reasoning 
applies to the fermions. Consider the action for any gaugino, 



d 5 x Re- 3Rky (-iXia n d n Xi) + / d 5 x Re~ iRky 



, d v -3/2Rk x 
-X 2 - i Xi + c.c. 

(51) 

The KK decomposition for these coupled differential equations diagonalizes 
the KK Dirac mass if 

(dy - 3/2Rk)X! = A 2 = 

Ai = (0 y - h/2Rk)X 2 = ^ 

This choice is again related to (JH?j) by ((7|). 
Brane Kinetic Terms 

We will eventually introduce a gauge invariant UV brane kinetic term for the 
SU(2)l gauge bosons to split the W ± and the charginos. Localized gauge 
kinetic terms modify the BCs, the definition of the scalar product and of the 
coupling constants. For example, such a term for the SU{2) L vectors 



S^S + njd'x 5(y)nR (-\ V m °V np *L a mn L a )j , (53) 
leads to the canonical normalization conditions for the W^, Z and photon 
JdyR(f L12 ( n \y) 2 + f R12{n) (y) 2 ) + nRKf Ll2 ^\0) 2 = 1 (54) 
f dy R (7 L3(n) (y) 2 + f m(n) (y) 2 + f x(n) (y) 2 ) + nRnf L3 ^ (0) 2 = 1 (55) 



2 / \ 2 N 
Oq \ ( do 



' R w +K) W + U 1 l=1 - (56) 



The Neumann BC d y A%(0) = becomes 

(d y + m\R 2 K) Aj(0) = 0. 
The dimensionless constant k is naturally ~ C(l)- 
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KK Decomposition 

The 5D fields are split into a coefficient bearing the 4D dependence, and the 
KK wavefunctions 

<P(x,y) = J2<P in) ^)f in \y)- (57) 

n 

The canonical normalization of the wavefunctions is J^dy Rf(y) 2 = 1 in flat 
space, but depends on the field in warped space because of the factor yfg 
and the different number of vielbeins e" For vectors, spinors and scalars 
respectively it is 

JdyRft{y) = l, jdyRe- 3Rk yfUy) = l, JdyRe^ f^(y) = 1 . (58) 

Any BC mixing two fields on the boundaries causes the two fields to belong 
to the same 4D KK tower, e.g. A^ n \x)f L{ - n \y) = A^ n \x)f R( - n \y) requires 
Ap (x) = Ap (x) to have the free action diagonal in the KK modes. 
Canonical normalization then means that the sum of the normalizations of 
all fields coupled in this manner should be unity. 
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